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ABSTRACT 

We use the Updated Zwicky Catalog of galaxies (Falco et al. 1999) to generate a 
catalog of groups, by means of a friend-of-friend algorithm. The groups identified have 4 
or more members and a surrounding density contrast, 5p/p = 80. They cover the region 
—4° < 5 < 90°, \b\ > 20° spanning over the radial velocity range of 2000 km s _1 < cz < 
15000 km s _1 . The total sample (hereafter GUZC) comprises 517 groups. The correlation 
length of the total sample is well fitted with a power law £(r) = (r/r ) 7 with parameters 
tq = 9.0 ± 0.4 h~ l Mpc and 7 = —1.67 ± 0.09 for values of r < 70 h~ 1 Mpc, declining to 
homogeneity at larger scales. 

Three subsamples defined by the range of group virial masses Ai were used to have 
their clustering properties examined throughout the autocorrelation function. We find 
an increase of the amplitude of the correlation function according to the group masses 
which extends the results of the ro — d c relation for galaxy systems at small d c . We obtain 
r = 9.5±0.5h~ 1 Mpc (7 = -1.81±0.12), r = 10.8±0.7h~ 1 Mpc (7 = -1.77±0.17), and 
r = U.l±1.2h~ 1 Mpc (7 = -1.65±0.22) for mass ranges 5 x 1O 12 M < M < 4xlO 14 M 
, 2 x 10 13 M o < M < 4 x 10 14 M Q and 5 x 1O 13 M < M < 4 x 1O 14 M respectively. 

For completeness we have also analyzed a sample of groups obtained from the Southern 
Sky Redshift Survey (da Costa et al.1998) in the range of virial masses 5 x 10 12 M Q < 
A4 < 4 x 1O 14 M to compare the results with those obtained from GUZC. The correlation 
function of this sample can be fitted with a power law with parameters r = 8.4 ± 
1.8 h^Mpc and 7 = -2.0 ± 0.7. 

Subject headings: cosmology: large-scale structure of universe — galaxies: clustering - 
galaxies: statistics 
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1. INTRODUCTION 



Th e two point correlation function, £(r), has 



been widely used to analyze the clustering prop- 
erties of galaxies and galaxy systems ( Peebles] 
1993) and references therein). Clusters and 



groups of galaxies can be used as suitable trac- 
ers of the large scale structure of the universe. 
In fact, several works have characterized the clus- 
tering properties of rich clusters of galaxies (see 
for instance |Croft et al. (199TJ| and |Abadi et al 



1998)|) concluding that the cluster autocorrela- 



tion function has a similar shape than the galaxy 
correlation function with a larger amplitude. The 
clustering properties of groups of galaxies were 
analyze d by |Jing & Zhang (1988]| for the CfA 
groups (|Geller fc Huchra (1983)| ), by |Maia fc da 
Costa 



1990 



et al. (1989 



for the S SRS group catalog ([Maia 
and by [Ramclla ct al. ( 1990 )| for 
groups of the CfA slices ( [Huchra et al. (199"0]D 
where a significant lower amplitude of £(r) is de- 
tected for these systems when compared to rich 
clusters, as expected in hierarchical scenarios of 
structure formation in the universe. However, the 
typical values found for the correlation length, r , 
for the groups examined by these authors are sig- 
nificantly different, perhaps due to the small size, 
and the different selection criteria of the sam- 



ples examined. In a more detailed work [lrasarti- 



IBattistoni et al. (1997) examine how the different 
criteria of assignment of galaxies into groups af- 
fect the determination of £(r). They claim that 
the differences found in previous works may be 
due to the distinct values of the density contrast 
adopted to generate the catalogs of groups, which 
resulted in groups with distinct physical proper- 
ties. 

Groups have lower richness and higher space 
density than the APM, Abell or X-ray selected 
clusters, allowing a detailed analysis of the be- 
havior of the tq at such regimes of space density. 
It worths to be mentioned that the relation be- 
tween r and the mean intergroup separation, d c , 
for values of d c lesser than 20h~ 1 Mpc has not been 
extensively explored. Moreover, the predictions 
of current hierarchical models differ from the uni- 
versal scaling law prop osed by [Bahcall and West 
(1992)| at large scales QAbadi et al. (1998)| and 
Croft et al. (1997)| ), as well as at the smaller 



scales where this work is focused. 

In this paper we analyze two new samples of 
groups of galaxies, the first one, taken from the 



Updated Zwicky Catalog (UZC) by [Falco et~aT 
(1999)| which contains redshift information for 



galaxies of the Catalog of Galaxies and Clusters 
of Galaxies ( Zwicky et al. (1961-68) ). The other 
group catalog, is derived from Southern Sky Red- 
shift Survey (SSRS2) by |da Costa et al. (1998) 



Both samples resulted, mainly, from the effort of 
the ON and CfA redshift surveys in the past 20 
years, but since UZC and SSRS2 come from dis- 
tinct galaxy catalogs, we treat both as indepen- 
dent samples, thus avoiding the mixture of each 
one biases, also allowing us, a cross-check of re- 
sults. 

Section 2 describes how the groups were defined 
and which are their general properties. In section 
3 we analyze the clustering properties of these 
groups using the autocorrelation function. A brief 
summary of our conclusions is presented in section 
4. 

2. GROUPS IN THE UPDATED ZWICKY CATALOG 
- GUZC AND SOUTHERN SKY REDSHIFT 
SURVEY - GSSRS2 



The recent release of the UZC of galaxies ( [Falco 
et al. (1999)|) , containing 19,369 objects with 



mzw < 15.5 and with 96% of completeness in 
redshift, and the SSRS2 ( |da Costa et al. (1D98)| ), 
containing 5369 galaxies with mg < 15.5 which 
covers a region of 1.70 sr of the southern celes- 
tial hemisphere. The large area of sky covered, 
the relative fainter magnitude limit, and the very 
good completeness in redshift of both catalogs, 
make them a powerful resource to perform sta- 
tistical analysis of the galaxy properties in the 
nearby Universe. We used the UZC and SSRS2 
to generate new catalogs of groups of galaxies to 
investigate their clustering properties by means of 
the correlation function. The algorithm adopted 
for the construction of the catalog of groups of 



galaxies is basically the one described by [Huchra 
|fc Geller (1982)| . Since we are interested in "real 
physical systems", the adopted friends-of-friends 
algorithm (described below), searches for possi- 
ble group member galaxies, keeping a fixed sur- 
rounding density contrast {Sp/p) relative to the 
mean density of galaxies in the catalog. A 5p/p = 
80 was adopted because it was demonstrated by 
Ramella et al. (1997)| to be the best compromise 



in identifying as many physical loose groups as 
possible and including all systems with high veloc- 
ity dispersion, but avoiding contamination of the 
catalog by pseudo-groups as well as groups with 
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interlopers (see als o discussion by [Nolthenius and 
Whh je (1987)1 and |Maia fc da Costa (1990)|) ~ 

A search for companions around galaxies is car- 
ried out taking into account projected separations 
satisfying 



2 sin 



7 12 



-Y- 



N 2 J H 

and with line-of-sight velocity differences, 



12 



v l -v 2 \<V L 



In the above expressions V = (Vi + V2V2, V\ and 
V2 are the radial velocities of the galaxies, and Qyi 
their angular separation. The quantities Dl and 
Vl are search parameters scaled according to the 
expressions below in order to take into account 
the variation in the sampling of the galaxy lumi- 
nosity function, <f)(M), with distance 



Di 



DnR 



Vl = V R 



where 



R 



A/12 



$(M)dM 



$(M)dM 



-1/3 



M Um = m Hm - 25 - 5 \og(V f /H 



M 



12 



m Hm -25-5\og[(V 1 + V 2 )/2H 



Dq is the selection parameter at a fixed fiducial 
radial velocity, Vf. Vl is scaled in the same way 
as Dl. 

To generate the UZC groups (GUZC), the 
adopted values for Do and Vo were 0.229/i~ 1 Mpc 
and 350 km s _1 respectively; the apparent mag- 
nitude limit, mu m = mz w = 15.5, and H Q = lOO/i 
km s _1 Mpc -1 . The groups present a surround- 
ing 5p/p = 80, for a |Schechter (1976)| luminos- 
ity function adopted for the total sample, param- 
eterized by 0* = 0.005 galaxies mag -1 Mpc -3 
M* 

(199-0)- 

The SSRS2 groups (GSSRS2) are obtained with 
the same Sp/p = 80, but for values for D and V 
0.352h~ 1 Mpc and 350 km s _1 respectively. The 
apparent magnitude, mu m = ms(o) = 15.5, while 
the luminosity function is parameterized by <fi* = 



l b(o) — -18.97 and a = —1.13 QRamella et al 



Mb® 



-19.40 



0.0137 galaxies mag 1 Mpc 
and a = -1.08 QMarzke et al. (1998)| ). 

We have also computed the virial masses of the 
groups combining the virial radius and the ve- 
locity dispersion M. cx R v cr 2 which provide suit- 
able estimates of group masses ( Hamella et al. 
(1997) ). According to the same authors, several 
groups with 3 or 4 members may be unreliable, 
although we should expect a low degree of con- 
tamination by pseudo-groups (as small as 10%) 
by adopting N > 4. 

The resulting group catalogs contain systems 
with at least 4 members and mean radial veloc- 
ities, V g < 15,000 km s- 1 (GUZC), and V g < 
12,000 km s _1 (GSSRS2) consisting in the largest 
sample of groups of galaxies homogeneously se- 
lected to the present. For analysis purposes, we 
removed groups in regions of galactic latitude 
\b\ < 20° and declination 5 > -4° for GUZC 
and for the GSSRS2 we only considered groups 
within the following limits: —40° < 5 < —2.5° 
and b < 40° (southern galactic cap) and 5 < 0° 
and b < +35° (northern galactic cap). This choice 
provides a sample free of strong obscuration from 
the Galaxy and also avoids possible lack of homo- 
geneity in the Falco catalog due to the not com- 
plete coverage of the Zwicky catalog in part of 
those regions. For both them, very rich groups 
containing more than 40 galaxies were also re- 
moved, to avoid the presence of rich clusters in 
our analysis. In addition to the above criteria, 
those groups with V g < 2,000 km s" 1 were also 
discarded from the analysis, to prevent any sig- 
nificant contribution of the group peculiar veloci- 
ties in the measured redshifts which were used to 
determine their distances. 

The projected distribution of the GUZC (panel 
a) and GSSRS2 (panel b) groups in equatorial 
coordinates is shown in figure 1. Histograms of 
the distributions of V g and group mass (M.) for 
the systems of the GUZC (panel a) and GSSRS2 
(panel b) are presented in figures 2 and 3 respec- 
tively. The velocity distribution of groups is sim- 
ilar to the one for galaxies, and a more detailed 
examination of the estimated masses for groups 
between 2000 km s _1 and 5000 km s" 1 confirms 
the reliability of this determination, in the sense 
that groups with lower masses are made up of few 
bright galaxies or, if they have a large number of 
members, they also contain a higher fraction of 
intrinsically less luminous objects. 
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3. ANALYSIS 



WPP . 



The lower richness and higher space densities 
of GUZC and GSSRS2 groups, compared to the 
APM or X-ray selected clusters allow the evalu- 
ation of the dependence of r with d c to be ex- 
amined for higher regimes of spatial densities in 
order to extend the available data of this relation 
and to compare with the r$ = 0.4 d c law ( |Bahcal] 
a nd West (19921) . 



The spatial correlation function is computed us- 
ing the estimator given by [Peebles (1993)| 



£(r) 



DD(r) N, 



R 



DR(r) N D 



1 



where DD(r) and DR(r) are the number of real 
and random pairs of groups at separation r and 
and Nr are the number of objects in the data 
and random catalogs, respectively. We use an 
equal weighting of each data/random point for 
simplicity given that error bars in our samples 
are significantly larger than the expected differ- 
ences between the estimates of £(r) from different 
weighting schemes according to the analysis by 
Ratcliffe et al. ( 1998 )|. The catalogs of random 



groups where generated within the same volume 
as the real groups, with a velocity distribution 
corresponding to the best third-degree polyno- 
mial fit to each V g distribution. The choice on 
the actual group distribution rather than that for 
the galaxies, is motivated by the fact that even 
though both present similar behavior, by using 
the group distribution it is possible to avoid bi- 
ases caused by the distribution of galaxies in low 
density regions devoid of richer groups. Errors 
were derived from the bootstrap resampling tech- 
nique developed by [Barrow et al. (1984)| with 30 
bootstrap samples. 

We estimated £ (r) for three samples, taking into 
account different intervals of group masses, A4, 
which represents, also, three distinct d c values. In 
order to deal with reasonable number of systems 
in the subsamples, the adopted mass limits are: 
5 x 1O 12 M < M < 4 x 1O 14 M , (397 groups, 
sample 1); 2 x 1O 13 M < M < 4 x 1O 14 M , 
(268 groups, sample 2); and 5 x 1O 13 M < M < 
4 x 1O 14 M , (167 groups, sample 3). The lower 



The estimated £(r) for the total sample of 
GUZC groups is shown if figure 4. The correlation 
function gives a positive signal up to 70h~ 1 Mpc 
and it has a similar shape to the one for the galax- 
ies. The results for the three constrained sub- 
samples are shown in figure 5, 6 and 7 for the 
range of masses corresponding to samples 1, 2 
and 3 respectively. By the inspection of figures 
4, 5, 6 and 7 it is possible to notice that the 
behavior of £(r) for all the four samples is con- 
sistent with a power-law. The amplitude of the 
correlation function presents a tendency to in- 
crease, as far as we have a predominance of more 
massive groups. We have used the Levenbergh- 
Marquardt method, which takes into account er- 
rors with minimum non-linear least-squares, to es- 
timate the best fitting parameters r$ and 7 in 
the power-law approximation (r/r ) 7 . We ob- 
tain r = 9.5 ± O.Sh^Mpc (7 = -1.81 ± 0.12), 
r = 10.8 ± 0.7h~ l Mpc (7 = -1.77 ± 0.17), and 
r = 14.1 ± \flhr x Mpc (7 = -1.65 ± 0.22) for 
subsamples 1, 2 and 3 respectively. 

In order to test the stability of the estimates 
of the correlation function against the existence 
of possible fake groups, we have removed from 
the total sample those groups with four mem- 
bers and mass exceeding 5 x 1O 13 M (48 groups). 
These massive groups with few members are sub- 
ject to suspicion and could strongly contribute to 
the expected 10% fraction of pseudo-groups as 
discussed by [Ramclla et al. (1997) . The corre- 
lation function for this restricted sample (r = 
(9.0 ± 0A)h~ 1 Mpc,-f = 1.72 ± 0.10) is indistin- 
guishable within errors from that of the total sam- 
ple, indicating that fake groups are not likely to 
generate strong systematic effects in our sample. 

Using the same procedure we have computed 
£(r) for the 104 groups in the GSSRS2 with 
the results displayed in figure 8. The best fit- 
ting parameters are r = 8.4 ± 1.8h _1 Mpc and 
7 = —2.0 ± 0.70. Due to the small number of 
objects in this sample we have not attempted di- 
visions into subsamples. 

A comparison of our results with those obtained 
by |Jing fc Zhang (1988)1 [Maia fc da Costa (19901 
IRamclla et al . (1990)1, and ITrasarti-Battistoni 



limit lin the first sample corresponds to the tvpi- et al. (1997)| is presented in Table 1. The r 



cal mass of the group with four L* members, the 
upper limit in the three samples was chosen to 
avoid groups with masses greater than the esti- 
mated mass of the Virgo cluster, ( [Ramella et al 



and 7 values they have found, are different when 
compared to the one we have found for the to- 
tal sample. The possible causes of the differences 
may be related to the small number of objects 



5 



in the previous samples, to the minimum num- 
ber of member galaxies assigned to groups, and 
also to the different Sp/p adopted for the defini- 
tion of groups, which results in the determination 
of groups with distinct physical properties. The 
different characteristic depths of the samples may 
also play an important role in the range of fit- 
ting £(r). In particular the samples analyzed by 
|Jing fc Zhang (1988) and by |Maia fc da Costa] 
( 199(0, which retain most of those characteristics 
described above, are the ones which present the 
most discrepant results. 

We have estimated the values of d c for the sub- 
samples in two different ways: (a) we compute 
the abundance n(> A4) of systems using |Bahcall| 
and ( Pen ( 1993^1 mass function, in the correspond- 
ing range of masses and estimate d c = t^ -1 / 3 ) . (b) 
we consider that the GUZC and GSSRS2 catalogs 
are bona-fide complete out to 5000 km s _1 and we 
fit a homogeneous distribution in the range 2000 
km s _1 < V^< 5000 km s _1 . Taken into account 



Ramerra~eL al. (1997)1 expect the loss of groups West (1992 



4. CONCLUSIONS 

We have analyzed the two point correlation 
function for two new catalogs of groups of galax- 
ies. The GUZC catalog is the largest homoge- 
neous sample to the present and provides a suit- 
able data set to analyze the clustering proper- 
ties of systems of galaxies of low richness. The 
group correlation function, for the total sample 
of groups, is well fitted to a power-law of the 
form £ = (r/r ) 7 , with 7 = —1.67 ± 0.09 and 
r = 9.0 ± 0.4 fr 1 Mpc for values of r < 70/r 1 
Mpc. A similar analysis for GSSRS2 catalog pro- 
duces r = 8.4 ± 1.8h -1 Mpc and 7 = -2.0 ± 0.7. 
These results differ from the findings of other au- 
thors mainly due to the different groups charac- 
teristics of the samples analyzed. 

The results for subsamples of groups with dif- 
ferent group mass intervals, show the tendency 
for higher values of ro at the small mean inter- 
group separations examined, when compared with 
the universal scaling law proposed by |Bahcall and 



at 50W km/s with at least 4 members less than 



20%, I so the sccund iiielhud ahuuld ulau piuvidcs 



a reliable estimate of d c . 

The resulting values for GUZC are: d c = 
[12] h^ 1 Mpc for the total sample, while d c = 
9 [14] h- 1 Mpc, d c = 15 [16] h' 1 Mpc and 
d c = 21 [22] h— 1 Mpc for subsamples 1, 2 and 

3 respectively. For the GSSRS2, d c = 9 [12] h' 1 
Mpc in the range of masses 5 x 1O 12 M < M. < 

4 x 10 14 M Q . The numbers in square brackets cor- 
respond to the estimates using method (b). Both 
methods give consistent values of d c , provided 
that the typical uncertainties in d c are ~ 2h~ l 



Mpcr- 

In figure 9 we plot the derived values of r as a 
function of d c for the different samples of groups 
considered. Also, are plotted in this figure the 



Correlation functions of APM clusters QCroft 
et al. (1997)|), X-ray confirm ed Abell clusters 
(XBACs, lAbadi et al. (1998)| ) and the results 
for groups presented in this paper strongly argue 
for a hierarchical scenario of structure formation 
in the universe such as cold dark matter models. 

This paper was partially supported by CON- 
ICET, SeCyT UNC and Fundacion Antorchas, 
Argentina. MAGM acknowledges CNPq grant 
301366/86-1. 

NOTE ADDED IN PROOF 
After this paper was submitted to ApJ, |Girardi 
et al. (2000)| examined £ (r) of groups relative to 



results for APM clusters according to |Croft et al 



(199f]]. All the group sets present r values higher 
than the universal scaling law tq = 0Ad c by Pah 



the one for galaxies using the UZC and SSRS2 
catalogs of galaxies. They find similar results to 
ours for a volume-limited sample of 139 groups 
(r = (8 ± l^Mpc; 7 = -1.9 ± 0.7). 



call and West (1992)1 in contrast to the clusters 



which have a lower amplitude than predicted by 
this law. Meanwhile, it is noticeable in this fig- 
ure the shape of the r — d c relation is similar to 
that expected in hierarchical models for structure 
formation such cold dark matter scenario dis- 
cussed by |Croft et al. (1997) . 
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Table 1 

Results of clustering analysis for several catalogs of groups. 



Sample 




r [h^Mpc] 


7 


d c [h^Mpc] 


# of groups 


Total Sample 




9.0 ±0.4 


-1.67 ±0.09 


- [12] 


517 


Sample 1 




9.5 ±0.5 


-1.81 ±0.12 


9 [14] 


397 


Sample 2 




10.8 ±0.7 


-1.77 ±0.17 


15 [16] 


268 


Sample 3 




14.1 ± 1.2 


-1.65 ±0.22 


21 [22] 


167 


GSSRS2 




8.4 ± 1.8 


-2.0 ±0.7 


9 [12] 


104 


Jing & Zhang (1988) 




3.5 


-1.77 




174 


Maia & da Costa (1990) 




3.2 


-1.82 


9.1 


87 


Ramella et al. (1990 




6.0 


-1.0 




128 


Trasarti-Battistoni et al. (1997) 


7.03 ± 0.6 


-1.21 ±0.04 




ps 200 



Quoted values of d c were calculated using Bahcall and Cen (1993) mass function, numbers in brackets 
correspond to the assumption of a homogeneous distribution in the range 2000 km s _1 



< V g < 5000 km s" 1 
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Fig. 1. — Projected distribution of GUZC (panel a) and GSSRS2 (panel b) groups. Solid lines correspond to the limit 
in galactic latitude |b| = 20° (GUZC) and |b| = 35° (limit of SSRS2 catalog) used in the statistical analysis. 
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Fig. 2. — Distribution of mean radial velocities (V g ) of GUZC (panel a) and GSSRS2 (panel b) groups. 
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Fig. 3. — Distribution of virial masses for GUZC (panel a) and GSSRS2 (panel b) groups. 
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Fig. 4. — Spatial autocorrelation function, £(r), for the total GUZC groups (517 objects). Error bars correspond 
to bootstrap resampling estimates of uncertainties. The solid line shows the power law fit (r/r ) 7 with r = 9.0 ± 0.4 
h^Mpc and 7 = -1.67 ±0.09. 
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Fig. 5.— The same as figure 4 for GUZC groups with 5 x 10 12 < M < 4 x 10 14 . 
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Fig. 6. — The same as figure 4 for GUZC groups with 2 x 10 13 < M < 4 x 10 
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Fig. 7.— The same as figure 4 for GUZC groups with 5 x 10 13 < M < 4 x 10 14 . 
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Fig. 8.— The same as figure 4 for GSSRS2 groups with 5 x 10 12 < M < 4 x 10 14 . 
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Fig. 9. — Correlation length, ro, as a function of the mean intergroup separation, d c . Circles correspond to the three 
GUZC subsamples of groups an alyzed, while triangles r epresents the values for the GSSRS2 sample. Open symbols 
correspond to d c estimates from Bahcall and Cen (1993) mass function, filled symbols to our estimates assuming com- 
pleteness in 2000 km s _1 < V g < 5000 km s _i . Squares show the ro — d c relation of APM clusters as described by Croft 
et al. 1997, while the line represents the universal scaling law r = 0Ad c . The error bars for our points in the figure 
correspond to those derived from the non linear least-squares fit. 



